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We model the spacetime of low-mass X-ray binaries with the Tomimatsu-Sato 
5 = 2 (TS2) metric and study the properties of the orbital and the epicyclic fre- 
quencies. The numerical analysis shows that the properties of the characteristic 
■ frequencies of oscillation do not differ qualitatively from those of the Kerr black 

hole. Estimates for the angular momenta of the three stellar mass black hole 
candidates GRO 1655-40, XTE 1550-564 and GRS 1915+105 are made with the 
application of the nonlinear resonance model. We find agreement between the 
5_i , predictions based on the 3 : 2 nonlinear resonance model for a TS2 background 



and the current estimates found in the literature. 
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1 Introduction 



Quasiperiodic oscillations (QPO) in the X-ray flux from low-mass (galactic) binaries 
has attracted considerable research interest recently due to their potential to be applied 
as a tool for probing the space-time around black holes and massive stars and testing 
gravity in strong field regime [U El |3_1 IH [5] . One of the major possible application of 
QPOs which makes them so attractive is the measurement of the angular momentum 
of the central object in X-ray binaries. There are very few methods for this measure- 
ment [6]: continuum fitting of thermal spectra [TJ, IHl |9] , determining the shape of the 
gravitationally redshifted wing of Fe line[T0] and QPOs [11]. Three types of QPO 
are usually observed: low-frequency QPOs (few tens of Hertz), intermediate frequency 
QPOs and high frequency QPOs (up to kilo Hertz). While the former two types of 
QPOs are believed to have astrophysical origin, i.e. they are related with the physics 
of the accretion disc, the high frequency QPOs are believed to depend on/reflect the 
structure of space-time in the vicinity of the central object - a black holes or a neutron 
star. Some properties of the QPOs support that idea: high-frequency QPOs depend 
very weakly on the X-ray flux and in several sources the frequencies of the QPOs occur 
in pairs whose ratio is 3 : 2. The physical mechanism of the production of these QPOs 
is not know. One of the major hypothesis is proposed by the nonlinear resonance 
model (NRM) according to which QPOs are related to the characteristic frequencies of 
oscillation - epicyclic and orbital - of test particles orbiting around the central object 
and the appearance of nonlinear resonances between them. This idea has initially been 
proposed by Aliev and Galtsov [H1[T5]. It has reappeared and further elaborated in the 
works of Abramowicz, Kluzniak and collaborators [161 El HE]- Simplicity makes the 
nonlinear resonance model very attractive, however, it has some significant deficiencies 
[20] . The model does not propose an excitation mechanism for the QPOs. The origin 
of the coupling between the different characteristic frequencies is not clear. The effect 
from astrophysical complication such as turbulence flow and magnetic instability have 
not been assessed and cannot be simply included in the model. All QPO models are 
purely dynamical and do not concern the emission mechanism of X-rays. A major dif- 
ficulty of the NRM is the dissonance between its predictions for the angular momenta 
of the observed black holes candidates and the measurements based on spectral contin- 
uum fitting [T9l |2"U] . Non of the different versions of the NRM can explain the observed 
angular momenta of the three black-hole candidates GRO 1655-40, XTE 1550-564 and 
GRS 1915+105. In the current paper we will address only the last one of the problems. 
According to Bambi [12] the dissonance between the predictions of the NRM and the 
measurements for the angular momenta mean that one of the following four possibili- 
ties must be true: i) The NRM is wrong ii) continuum-fitting method does not provide 
reliable estimates of the angular momenta Hi) both techniques do not work correctly 
iv) the central objects of the observed stellar-mass black hole candidates are not Kerr 
black holes. Here we will consider only the fourth possibility. 

The effect of deviations from Kerr space-time on the QPOs in the frame of RNM 
has already been studied in several papers [13] . [T2] , [5], [2TJ. In all of these cases the 
deviation from Kerr space-time is described by an additional continuous free parameter 
(e, tidal charge, etc.) which leads to degeneracy of the results in the sense that for fixed 
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of the central object mass and value of the deviation parameter there is a whole interval 
of values of the angular momentum which are admissible by the observations. Due to 
this degeneracy, in order to be able to determine the value of the deviation parameter 
we need independent measurements of the mass and the angular momentum. One of 
the few exact solutions that describe the exterior field of a stationary rotating axially 
symmetric objects is the Tomimatsu-Sato solution. The deviation parameter 5 in this 
solution takes discrete values (5=1 corresponding to Kerr solution 5 = 2 corresponding 
to an exotic object) which allows the degeneracy to be avoided . The aim of the present 
paper is to apply the NRM to the space-time of Tomimatsu-Sato with 5 = 2 (TS2), 
to study the properties of the characteristic frequencies and to check if the dissonance 
problem can be resolved in such background. The paper is organized as follows. The 
Tomimatsu-Sato spacetime is briefly presented in Section [2j The formulas that we use 
for the calculation of the epicyclic frequencies are obtained in Appendix A and applied 
in Section [3j In Section H] NRM model is sketched. The results for the estimates of the 
angular momenta of the three black hole candidates are in Section [5j The sum-up of 
the results is in Section [6] 



2 Tomimatsu-Sato space-time 

The Kerr- Tomimatsu-Sato family is usually presented in the canonical Weyl-Papapetrou 
form of the metric for stationary, axisymmetric space-time [22] 

ds 2 = -f(dt - ud<f)) 2 + f-^e^idp 2 + dz 2 ) + p 2 d(j)\ (1) 

where [23] 

f = t - = ^d-y 2 )c, ^ = -_d_ . (2) 



The functions A, B and C are polynomials of the prolate spheroidal coordinates x and 

y defined by 

p = cr a/ (x 2 - 1)(1 - y 2 ), z = axy. (3) 

In explicit form for the Kerr space-time 5 = 1 these polynomials are 

A = p 2 (x 2 -l)-q 2 (l-y 2 ), B = (px + l) 2 + q 2 y 2 , C = -px - 1. (4) 
In the Tomimatsu-Sato case 5 = 2 , A, B and C are given by 

A = p\x 2 - If + q\l - y 2 f 

-2p 2 q 2 (x 2 - 1)(1 - y 2 ){2(x 2 - l) 2 + 2(1 - y 2 ) 2 + 3(x 2 - 1)(1 - y 2 )} 
B = {p 2 (x 2 + l)(a; 2 - 1) - q 2 (y 2 + 1)(1 - y 2 ) + 2px(x 2 - l)} 2 

+Aq 2 y 2 {px{x 2 - 1) + (px + 1)(1 - y 2 )} 2 (5) 
C = -p 3 x(x 2 - l){2(x 2 + l)(x 2 - 1) + (x 2 + 3)(1 - y 2 )} 

-p 2 [x 2 - l){Ax 2 (x 2 - 1) + (3x 2 + 1)(1 - y 2 )} + q 2 (px + 1)(1 - y 2 f. 
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The Kerr-TS family has three free parameters a, p and q which are related to the mass 
M and the angular momentum J of the central object as 

p 2 + q 2 = l, ° = ^f, J = M\ (6) 

The properties of the TS2 space-time have been studied in several papers [24"ll28| |29~ 1 130] . 
A simple sketch of the structure of TS2 can be found in [31]. We will just briefly 
mention its main properties. The TS2 space-time has two spherically symmetric Killing 
horizons symmetrically distributed with respect to the equatorial plane at p = and 
z = a (x = 1 and y = 1) connected by a rod-like conical singularity which carries 
all the gravitational mass. The Killing horizons are everywhere regular except at 
the points of connection with the rod-like singularity. In the equatorial plane there 
is a ring singularity at the roots of B(x, y = 0) = with zero Komar mass and 
infinite circumference. The metric component diverge on the ring singularity and 
is negative between the ring singularity and the central rod-like singularity. The latter 
means that there is causality violation and closed time-like curves in that region. 

The formula for the quadrupole momentum of the Kerr-Tomimatsu-Sato family has 
been correctly derived in [30] 

Q = -M 3 (^V + g 2 ) • (7) 

Due to its peculiar causal structure the TS2 space-time is discarded as a candidate 
to describe the final stage of gravitational collapse. The more appealing physical in- 
terpretation is the one given by Gibbons and Russell-Clark [23] that it describes the 
exterior field of a specific star whose prolateness is bigger than that of the Kerr black 
hole. The stability of TS2 has not been studied up to now. 



3 Characteristic epicyclic and orbital frequencies 

The epicyclic frequencies of relativistic compact objects to our knowledge have been 
derived for the first time in [14J. General formulas for the vertical and radial epicyclic 
frequencies can be found in [25] and [26J. In appendix |A] of the current paper we present 
a slightly different derivation. In this derivation the decoupling of the vertical and the 
radial modes is not assumed but rather follows from the reflectional symmetry of the 
background metric (For more details on reflectional symmetry see [27]-)- The vertical 
and radial epicyclic frequencies, respectively, are v 2 p = V p p /{2n) and v\ = V z z /(2tt). 
The matrix T> is defined in formulas (1291) and (I30I) form Appendix [A] For the particular 
background metric, the TS2 spacetime, the off-diagonal components vanish, T> z p = 
and T> p z = 0. For the Keplerian frequency Qk we apply the formula 

Q _ d ^> _ ~ d P 9t <f> ± \l ( d p9*<i>) 2 ~ ( 9 p9tt) {d P 9w) , g . 

Plus "+" refers to prograde orbits and minus "— ", to retrograde. For the calculation of 
the characteristic frequencies of the TS2 space-time we apply a package linear algebra 
in Maple. The formulas are too cumbersome to be displayed explicitly. 
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In some of the particular cases in which the resonance model was applied that we 
cited in the Introduction it was found that the characteristic frequencies may have 
rather odd behavior qualitatively different from the case of Kerr black hole. Vertically 
unstable modes have been reported in [2TJ [12], multiple (rather than one) extrema of 
the radial frequencies - in [321 |2T) . violation of the Uf- > i/g > v r inequality, where u$ 
and v T are the frequencies of the vertical and the radial modes respectively, has been 
shown to exist in [121 E21 EI], absence of radially unstable modes has been found in 

m 

Unlike the cases discussed above there are no peculiarities in the properties of the 
epicyclic frequencies of the TS2 space-time. Qualitatively their behavior is the same 
as in the case of Kerr black hole. The radial frequency v p has only one extremum 
(maximum). It increases as the singularity is approached, passes through a maximum 
and becomes zero at p = p/sccH before the singularity is reached. The frequency v p 
is presented on the left panel of figure (CQ) for q = — 0.9; 0; 0.9. The same behavior 
is observed for the whole interval of admissible values of q. As we can see from the 
representative cases shown on the figure for positive angular momenta of the central 
object the frequencies are higher than for the negative angular momenta. The verti- 
cal epicyclic frequency v z and the orbital frequency Uk increase monotonically as the 
singularity is approached (See the right panel of figure (jTJ and figure (J2J) ). For these 
two frequencies the effect of q is opposite in comparison to the radial frequency case. 
Negative angular momenta invoke higher frequencies while positive angular momenta 
- lower frequencies. All three frequencies u p , v z and vk are displayed on figure ([3]) for 
q = 0.9. For corotating orbits v p < v z < vr. As we have already mentioned the TS2 
space-time is a significant deformation of Kerr space-time. The properties of these two 
space-times differ significantly even when q = 0. On figure (jlj) v p for q = is shown and 
as it can be seen v p is lower in the TS2 case. The radius of the TS2 ISCO is higher than 
the one of Kerr. The differences of v z and vk in the two cases are indistinguishable on 
the graphics. 

4 Nonlinear resonance model for high frequency 
QPOs 

According to the NRM resonances occur at the orbits for which the relation 

mu p (q,r) = nv z {q,r) (9) 

is satisfied for integer "m" and "n". The values of "m" and "n" vary in the different 
resonance models. When "g" is fixed ([9j) can be solved for the radii of the resonance 
orbits r nm . For the observed twin peak high-frequency QPOs the upper vy and the 
lower ul frequency are in 3 : 2 ratio so within the model they are expressed in the 
following way 

Vu = miVp + m 2 v z , = n x v p + n 2 v z , (10) 

1 ISCO stands for innermost stable circular orbit. 
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p p 

Figure 1: The radial and the vertical epicyclic frequencies, v p and u z , for q = 
-0.9; 0; 0.9. 
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Figure 3: The radial and the vertical epicyclic frequencies v p and v z and the orbital 
frequency u K for q = 0.9. 




— TS2 

— Kerr 



Figure 4: A comparison between the Kerr and TS2 radial frequencies v p for q = 0. 
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Model 


GRO 1655-40 


XTE 1550-564 


GRS 1915+105 


M/Mq 


6.30 ±0.27 


9.1 ±0.6 


14.0 + 4.4 


vu, [Hz] 


450 ±3 


276 ±3 


168 + 3 


vl, [Hz] 


300 ±5 


184 ±5 


113 + 5 



Table 1: Stellar- mass BH candidates in microquasars with a measurement of the mass 
and two observed high-frequency QPOs [321 EH 133] . 

where the epicyclic frequencies are evaluated at the resonance orbits. The small integers 
{mi, m,2, ni, 712} are chosen appropriately so the vu : vl = 3 : 2 ratio is satisfied. For 
the forced oscillation model (777,77) = (3, 2). In the parametric resonance model either 
(m, n) = (3, 1) or (777, n) = (2, 1). 

In the keplerian resonance model the resonance condition is, respectively, 

mu p (q,r) = nv K {q,r), (11) 

where (777, n) = (3, 2); (3, 1); (2, 1). The upper and the lower twin peak QPO frequencies 
are 

vu — m\v p + m^v K , u L — niu p + n2VK- (12) 

Again the radial epicyclic frequency and the orbital frequency are evaluated at the 
resonance orbits. 

5 Estimates for the angular momenta of black-hole 
candidates in microquasars 

In this section we apply different nonlinear resonance models to estimate the angular 
momenta of the central object in three microquasars GRO 1655-40, XTE 1550-564 and 
GRS 1915+105. For these black holes candidates we have also dynamical measurement 
of the mass. For all three of them two high-frequency QPOs have been observed and 
the ratio between the upper vy and the lower frequency Vl is 3 : 2 (See Table. [[]). 
Two high-frequency QPOs have been observed also in the X-ray binary H 1743-322, 
however, there is no measurement for its mass. 

The results from the numerical calculation for the estimates of the angular momenta 
of the microquasars based on the models in which nonlinear resonance occurs between 
the radial and vertical epicyclic frequencies are presented in Table. [2] and for resonance 
between the radial epicyclic frequency and the orbital frequancy - in Table. [3J 
In general we can say that when the central object of the microquasar is modeled 
by the central object of TS2 space-time the estimates for the angular momenta are 
lower than in the case when the central object is modeled by a Kerr black hole. The 
estimates in the latter case can be found in [17] . One of the most significant problems 
of the nonlinear resonance model is the discrepancy between the estimates for the 
angular momenta based on the observed QPOs and the estimates obtained thought 
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Model GRO 1655-40 XTE 1550-564 GRS 1915+105 



3:2 0.63 - 


- 0.68 0.41 - 


- 0.55 -0.11 - 


- 0.68 


2:1 0.40 - 


- 0.49 0.21 - 


- 0.38 -0.40 - 


- 0.55 


3:1 0.58 - 


- 0.66 0.42 - 


- 0.57 -0.12 - 


- 0.71 



Table 2: Angular momentum estimates for the three microquasars GRO 1655-40, 
XTE 1550-564 and GRS 1915+105 from the models in which the resonance occurs 
between the radial and vertical epicyclic frequencies. 



Model GRO 1655-40 XTE 1550-564 GRS 1915+105 



3:2 0.62 - 


- 0.71 0.45 - 


- 0.61 -0.14 - 


- 0.76 


2:1 0.40 - 


- 0.49 0.21 - 


- 0.38 -0.40 - 


- 0.55 


3:1 0.54 - 


- 0.61 0.39 - 


- 0.52 -0.10 - 


- 0.65 



Table 3: Angular momentum estimates for the three microquasars GRO 1655-40, XTE 
1550-564 and GRS 1915+105 from the models in which the resonance occurs between 
the radial epicyclic frequency and the orbital (Keplerian) frequency. 

analysis of the observed thermal spectra of the objects' accretion disks. In non of 
the resonance models there is agreement between QPOs and thermal spectra for all 
three microquasars. For example, the measured value of angular momentum of GRS 
1915+105 coming from the continuum-fitting method is in agreement with the 3 : 2 
epicyclic resonance model but the value predicted by the same resonance model for 
GRO 1655-40 does not match the spectral fitting estimate. It would be interesting to 
know whether such discrepancy is still present when TS2 solution is used to describe 
the space-time in the vicinity of the microquasar. In order to be able to legitimately 
compare between the estimates of the resonance model of QPOs and those of the 
thermal spectra fitting the experimental data in the later case should be reanalyzed 
with TS2 background. In the absence of such analysis we could compare the estimates 
of the resonance model for TS2 with the estimates of the spectral continuum-fitting 
and other methods for Kerr. The continuumm-fitting method [8] gives the following 
values for GRO 1655-40, a G (0.65 + 0.75). The angular momentum of XTE 1550-564 
has been estimated in j9] to be a G (0.29 + 0.52). These values have been obtained 
through the combination of the results from the continuum-fitting method and the Fe- 
line method. In literature two conflicting estimates for GRS 1915+105 can be found 
(See also |36j). The first one is given in [7] - a > 0.98. The estimate given in [37] is 
a ~ 0.7. From tables [2] and [3] we can see that the 3 : 2 Keplerian resonance is consistent 
with these values for all three microquasars if the alternative estimate for the angular 
momentum of GRS 1915+105 [37] is chosen. 
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6 Sum-up 



In the current paper we have studied the properties of the three characteristic frequen- 
cies of a particle propagating alone a circular orbit in the in the equatorial plane of 
Tomimatsu-Sato space-time - the radial u p and the vertical v z epicyclic frequencies and 
the orbital (Keplerian) frequency vk ■ For the calculation of the frequencies we have fol- 
lowed a variation of the method proposed [2|. The formulas are presented in appendix 
A of the current paper. The numerical analysis shows that for the whole interval of 
admissible values of the free parameter q corresponding to angular momentum per unit 
mass the quantitative behavior of the frequencies are the same as in the case of Kerr 
black hole. We have also applied the resonance mode to estimate the angular momenta 
of three microquasars GRO 1655-40, XTE 1550-564 and GRS 1915+105. The central 
object here is modeled by the central object in the Tomimatsu-Sato 5 = 2 space-time. 
The results are presented in Table. [2] and Table. [31 The estimated angular mementa of 
the three microquasars are considerably lower than in the case when the central object 
is modeled by a Kerr black hole. 

The 3 : 2 keplerian NRM seems to be in agreement with the some of the current 
estimates for the angular momenta of three black hole candidates. However, future 
studies should take into account that the continuum-fitting method is model dependent 
and in the cases cited here it was supposed that the central object is a Kerr black hole. 

In conclusion, the results presented here do not exclude the Tomimatsu-Sato 5 = 2 
space-time as a model of the space-time of microquasars. 

A Derivation of the formula for the epicyclic fre- 
quencies 

x^ + T^x 9 ^ = (13) 

Sx" + 2Y t " pa x p 5x a + a x a xHx° = (14) 

We are interested in circular orbits in space-times with two Killing vectors and ^ 
- stationary (or static in particular) and axially symmetric . Due to the invariance 
of the space-time with respect to translations in time and rotations with respect to 
an axis the law of conservation of energy and the law of conservation of angular mo- 
mentum hold for point particles. We will use adapted coordinates in which the two 
Killing vectors have the form £ f = d/dt and ^ = d/d<p. The solutions whose QPOs 
we will study are presented Weyl-Papapetrou (t, p, z, 0) or Boyer-Lindquist (t, r, 9, 0) 
coordinates. The derivation of the epicyclic frequencies presented here is valid also if 
prolate spheroidal coordinates are used to present an axially symmetric solution. By 
flat orbits here we mean orbits that remain in a plane and in particular in the equa- 
torial plane. In the coordinate systems we have chosen these orbits are described as 
p = const and z = const or alternatively r = const and 9 = const. For convenience 
we will adopt the following notation. Small Latin letters a,b,c,d,s,m,n.... will be 
used to denote the r and 9 (respectively the p and z) coordinates. Capital Latin let- 
ters A, B, C, D, S, M, N.... represent the other two coordinates, those related with the 
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energy and angular momentum conservation laws, t and (j) respectively. We will also 
term them cyclic coordinates since the lagrangian of point particles does not depend on 
them. Greek letter will denote all coordinates a, (3, 7, 5, a, /1, v... = t, r, 9, (t, p, z, 0). 
In that notation flat, circular orbits are expressed as x a = where the dot denotes the 
derivative with respect to the affine parameter of the orbits - A. Due to stationarity 
and axial symmetry for the metric we have g M „ = g^x" 1 ). It can be proved that all 
quantities that are function of the x a coordinates only remain constant on the chosen 
type of orbits. Let / = f(x a ). Then 



since x a = 0. In particular 



/ = £ = f,.i° = 0, (15) 

(16) 

For the particular case of flat, circular orbits 'x a = while in general for the cyclic 
coordinates 'x A = holds. So 

x a = 0. (17) 
Due to (fTS]) and ([IT]) . dHJ) can be rewritten as 



dX 

Taking into account that 



' ' (Sx" + 2T^^5x a ) = -T\ Pta x a x p 5x° 



r% = \^ (d a9$<7 + dpg w - d a g afS ) (19) 



we can easily show that some of the Christoffel symbols vanish. 



f m ab = -\g Ms d s 9AB = 0, r m AB = - l -g ms d s9AB ± o ; 

r\ = \g MP d b g A P ^ 0, T\ b = l -g mP d b g AP = 0. (20) 

We have taken into account that g A B 7^ for any A and B, g ab 7^ for a = b, 
9Ab = = g a B and for the inverse metric g AB 7^ for A any B, g ab 7^ for a = b, 
9 Ab = = g aB . 

The M components of (I14p simplifies considerably. From (JTSj) and (|20|) 

^ (5x M + 2T M Bs x B 5x s ) =0 (21) 

so 

5x M + 2T M Bs x B 5x s = const = 0. (22) 

Setting the constant to zero is equivalent to choosing the difference between initial 
phases of 5x M and 5x a to be it/2. Since (d/dt) = ^(d/dX) from (122]) it follows that 

^- + 2 (rl + r^ K ) 5x s = 0. (23) 
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f^K is the Keplerian frequency 



= t = d 4 = i (24) 
t dt E K ' 



L and E are the angular momentum and the energy of the particle, respectively. Now 
let us obtain the equation that governs the m component of the perturbations 



5x m + 2Y m BP x B 5x p + T m AB s x A x B 5x s = 0. (25) 



It contains both the Sx m and the Sx perturbations. The Sx can be decoupled if ( 1221) 
is used 

5x m + 2T m BP x B (-2T p As x A Sx s ) + T m AB s x A x B 5x a = 0. (26) 

Finally, 

5x m + F m ABs x A x B 5x s = 0, (27) 
where we have introduced the following matrix 

F m ABs = T m ABtS -AT m AP T p Bs . (28) 
Dividing the upper equation in t we obtain 

+ {F m tts + 2F™ s n K + F^Ql) 5x s = 0, (29) 
d 2 5x m 



dt 2 tts 

or 

(lf2 +V m s 5x s = 0. (30) 

In the general case the horizontal ( r or p) and the vertical (8 or z) perturbations are 
coupled. When the "elasticity" matrix T> is diagonal the epicyclic frequencies are 

& m = V m ml m = r 1 9 (p,z). (31) 

No summation over the repeated index m is implied. 
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